Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows.

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.

A pan-influenza antibody inhibiting neuraminidase via receptor mimicry.

Rapidly evolving influenza A viruses (IAVs) and influenza B viruses (IBVs) are major causes of recurrent lower respiratory tract infections. Current influenza vaccines elicit antibodies predominantly to the highly variable head region of haemagglutinin and their effectiveness is limited by viral drift1 and suboptimal immune responses2. Here we describe a neuraminidase-targeting monoclonal antibody, FNI9, that potently inhibits the enzymatic activity of all group 1 and group 2 IAVs, as well as Victoria/2/87-like, Yamagata/16/88-like and ancestral IBVs. FNI9 broadly neutralizes seasonal IAVs and IBVs, including the immune-evading H3N2 strains bearing an N-glycan at position 245, and shows synergistic activity when combined with anti-haemagglutinin stem-directed antibodies. Structural analysis reveals that D107 in the FNI9 heavy chain complementarity-determinant region 3 mimics the interaction of the sialic acid carboxyl group with the three highly conserved arginine residues (R118, R292 and R371) of the neuraminidase catalytic site. FNI9 demonstrates potent prophylactic activity against lethal IAV and IBV infections in mice. The unprecedented breadth and potency of the FNI9 monoclonal antibody supports its development for the prevention of influenza illness by seasonal and pandemic viruses.

Towards increased accuracy and reproducibility in SARS-CoV-2 next generation sequence analysis for public health surveillance.

During the COVID-19 pandemic, SARS-CoV-2 surveillance efforts integrated genome sequencing of clinical samples to identify emergent viral variants and to support rapid experimental examination of genome-informed vaccine and therapeutic designs. Given the broad range of methods applied to generate new viral genomes, it is critical that consensus and variant calling tools yield consistent results across disparate pipelines. Here we examine the impact of sequencing technologies (Illumina and Oxford Nanopore) and 7 different downstream bioinformatic protocols on SARS-CoV-2 variant calling as part of the NIH Accelerating COVID-19 Therapeutic Interventions and Vaccines (ACTIV) Tracking Resistance and Coronavirus Evolution (TRACE) initiative, a public-private partnership established to address the COVID-19 outbreak. Our results indicate that bioinformatic workflows can yield consensus genomes with different single nucleotide polymorphisms, insertions, and/or deletions even when using the same raw sequence input datasets. We introduce the use of a specific suite of parameters and protocols that greatly improves the agreement among pipelines developed by diverse organizations. Such consistency among bioinformatic pipelines is fundamental to SARS-CoV-2 and future pathogen surveillance efforts. The application of analysis standards is necessary to more accurately document phylogenomic trends and support data-driven public health responses.

Imprinted antibody responses against SARS-CoV-2 Omicron sublineages.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron sublineages carry distinct spike mutations resulting in escape from antibodies induced by previous infection or vaccination. We show that hybrid immunity or vaccine boosters elicit plasma-neutralizing antibodies against Omicron BA.1, BA.2, BA.2.12.1, and BA.4/5, and that breakthrough infections, but not vaccination alone, induce neutralizing antibodies in the nasal mucosa. Consistent with immunological imprinting, most antibodies derived from memory B cells or plasma cells of Omicron breakthrough cases cross-react with the Wuhan-Hu-1, BA.1, BA.2, and BA.4/5 receptor-binding domains, whereas Omicron primary infections elicit B cells of narrow specificity up to 6 months after infection. Although most clinical antibodies have reduced neutralization of Omicron, we identified an ultrapotent pan-variant-neutralizing antibody that is a strong candidate for clinical development.

Antibody therapy reverses biological signatures of COVID-19 progression.

Understanding who is at risk of progression to severe coronavirus disease 2019 (COVID-19) is key to clinical decision making and effective treatment. We study correlates of disease severity in the COMET-ICE clinical trial that randomized 1:1 to placebo or to sotrovimab, a monoclonal antibody for the treatment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (ClinicalTrials.gov04545060). Laboratory parameters identify study participants at greater risk of severe disease, including a high neutrophil-to-lymphocyte ratio (NLR), a negative SARS-CoV-2 serologic test, and whole-blood transcriptome profiles. Sotrovimab treatment is associated with normalization of NLR and the transcriptomic profile and with a decrease of viral RNA in nasopharyngeal samples. Transcriptomics provides the most sensitive detection of participants who would go on to be hospitalized or die. To facilitate timely measurement, we identify a 10-gene signature with similar predictive accuracy. We identify markers of risk for disease progression and demonstrate that normalization of these parameters occurs with antibody treatment of established infection.

Predicting the mutational drivers of future SARS-CoV-2 variants of concern.

SARS-CoV-2 evolution threatens vaccine- and natural infection-derived immunity as well as the efficacy of therapeutic antibodies. To improve public health preparedness, we sought to predict which existing amino acid mutations in SARS-CoV-2 might contribute to future variants of concern. We tested the predictive value of features comprising epidemiology, evolution, immunology, and neural network-based protein sequence modeling, and identified primary biological drivers of SARS-CoV-2 intra-pandemic evolution. We found evidence that ACE2-mediated transmissibility and resistance to population-level host immunity has waxed and waned as a primary driver of SARS-CoV-2 evolution over time. We retroactively identified with high accuracy (area under the receiver operator characteristic curve, AUROC=0.92-0.97) mutations that will spread, at up to four months in advance, across different phases of the pandemic. The behavior of the model was consistent with a plausible causal structure wherein epidemiological covariates combine the effects of diverse and shifting drivers of viral fitness. We applied our model to forecast mutations that will spread in the future and characterize how these mutations affect the binding of therapeutic antibodies. These findings demonstrate that it is possible to forecast the driver mutations that could appear in emerging SARS-CoV-2 variants of concern. We validate this result against Omicron, showing elevated predictive scores for its component mutations prior to emergence, and rapid score increase across daily forecasts during emergence. This modeling approach may be applied to any rapidly evolving pathogens with sufficiently dense genomic surveillance data, such as influenza, and unknown future pandemic viruses.

Broadly neutralizing antibodies overcome SARS-CoV-2 Omicron antigenic shift.

The recently emerged SARS-CoV-2 Omicron variant encodes 37 amino acid substitutions in the spike protein, 15 of which are in the receptor-binding domain (RBD), thereby raising concerns about the effectiveness of available vaccines and antibody-based therapeutics. Here we show that the Omicron RBD binds to human ACE2 with enhanced affinity, relative to the Wuhan-Hu-1 RBD, and binds to mouse ACE2. Marked reductions in neutralizing activity were observed against Omicron compared to the ancestral pseudovirus in plasma from convalescent individuals and from individuals who had been vaccinated against SARS-CoV-2, but this loss was less pronounced after a third dose of vaccine. Most monoclonal antibodies that are directed against the receptor-binding motif lost in vitro neutralizing activity against Omicron, with only 3 out of 29 monoclonal antibodies retaining unaltered potency, including the ACE2-mimicking S2K146 antibody1. Furthermore, a fraction of broadly neutralizing sarbecovirus monoclonal antibodies neutralized Omicron through recognition of antigenic sites outside the receptor-binding motif, including sotrovimab2, S2X2593 and S2H974. The magnitude of Omicron-mediated immune evasion marks a major antigenic shift in SARS-CoV-2. Broadly neutralizing monoclonal antibodies that recognize RBD epitopes that are conserved among SARS-CoV-2 variants and other sarbecoviruses may prove key to controlling the ongoing pandemic and future zoonotic spillovers.

Broad betacoronavirus neutralization by a stem helix-specific human antibody.

The spillovers of betacoronaviruses in humans and the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants highlight the need for broad coronavirus countermeasures. We describe five monoclonal antibodies (mAbs) cross-reacting with the stem helix of multiple betacoronavirus spike glycoproteins isolated from COVID-19 convalescent individuals. Using structural and functional studies, we show that the mAb with the greatest breadth (S2P6) neutralizes pseudotyped viruses from three different subgenera through the inhibition of membrane fusion, and we delineate the molecular basis for its cross-reactivity. S2P6 reduces viral burden in hamsters challenged with SARS-CoV-2 through viral neutralization and Fc-mediated effector functions. Stem helix antibodies are rare, oftentimes of narrow specificity, and can acquire neutralization breadth through somatic mutations. These data provide a framework for structure-guided design of pan-betacoronavirus vaccines eliciting broad protection.

Broad sarbecovirus neutralization by a human monoclonal antibody.

The recent emergence of SARS-CoV-2 variants of concern1-10 and the recurrent spillovers of coronaviruses11,12 into the human population highlight the need for broadly neutralizing antibodies that are not affected by the ongoing antigenic drift and that can prevent or treat future zoonotic infections. Here we describe a human monoclonal antibody designated S2X259, which recognizes a highly conserved cryptic epitope of the receptor-binding domain and cross-reacts with spikes from all clades of sarbecovirus. S2X259 broadly neutralizes spike-mediated cell entry of SARS-CoV-2, including variants of concern (B.1.1.7, B.1.351, P.1, and B.1.427/B.1.429), as well as a wide spectrum of human and potentially zoonotic sarbecoviruses through inhibition of angiotensin-converting enzyme 2 (ACE2) binding to the receptor-binding domain. Furthermore, deep-mutational scanning and in vitro escape selection experiments demonstrate that S2X259 possesses an escape profile that is limited to a single substitution, G504D. We show that prophylactic and therapeutic administration of S2X259 protects Syrian hamsters (Mesocricetus auratus) against challenge with the prototypic SARS-CoV-2 and the B.1.351 variant of concern, which suggests that this monoclonal antibody is a promising candidate for the prevention and treatment of emergent variants and zoonotic infections. Our data reveal a key antigenic site that is targeted by broadly neutralizing antibodies and will guide the design of vaccines that are effective against all sarbecoviruses.

SARS-CoV-2 immune evasion by the B.1.427/B.1.429 variant of concern.

A novel variant of concern (VOC) named CAL.20C (B.1.427/B.1.429), which was originally detected in California, carries spike glycoprotein mutations S13I in the signal peptide, W152C in the N-terminal domain (NTD), and L452R in the receptor-binding domain (RBD). Plasma from individuals vaccinated with a Wuhan-1 isolate-based messenger RNA vaccine or from convalescent individuals exhibited neutralizing titers that were reduced 2- to 3.5-fold against the B.1.427/B.1.429 variant relative to wild-type pseudoviruses. The L452R mutation reduced neutralizing activity in 14 of 34 RBD-specific monoclonal antibodies (mAbs). The S13I and W152C mutations resulted in total loss of neutralization for 10 of 10 NTD-specific mAbs because the NTD antigenic supersite was remodeled by a shift of the signal peptide cleavage site and the formation of a new disulfide bond, as revealed by mass spectrometry and structural studies.

SARS-CoV-2 RBD antibodies that maximize breadth and resistance to escape.

An ideal therapeutic anti-SARS-CoV-2 antibody would resist viral escape1-3, have activity against diverse sarbecoviruses4-7, and be highly protective through viral neutralization8-11 and effector functions12,13. Understanding how these properties relate to each other and vary across epitopes would aid the development of therapeutic antibodies and guide vaccine design. Here we comprehensively characterize escape, breadth and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD). Despite a trade-off between in vitro neutralization potency and breadth of sarbecovirus binding, we identify neutralizing antibodies with exceptional sarbecovirus breadth and a corresponding resistance to SARS-CoV-2 escape. One of these antibodies, S2H97, binds with high affinity across all sarbecovirus clades to a cryptic epitope and prophylactically protects hamsters from viral challenge. Antibodies that target the angiotensin-converting enzyme 2 (ACE2) receptor-binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency. Nevertheless, we also characterize a potent RBM antibody (S2E128) with breadth across sarbecoviruses related to SARS-CoV-2 and a high barrier to viral escape. These data highlight principles underlying variation in escape, breadth and potency among antibodies that target the RBD, and identify epitopes and features to prioritize for therapeutic development against the current and potential future pandemics.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants.

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

Transfer transcriptomic signatures for infectious diseases.

The modulation of the transcriptome is among the earliest responses to infection. However, defining the transcriptomic signatures of disease is challenging because logistic, technical, and cost factors limit the size and representativeness of samples in clinical studies. These limitations lead to a poor performance of signatures when applied to new datasets. Although the study focuses on infection, the central hypothesis of the work is the generalization of sets of signatures across diseases. We use a machine learning approach to identify common elements in datasets and then test empirically whether they are informative about a second dataset from a disease or process distinct from the original dataset. We identify sets of genes, which we name transfer signatures, that are predictive across diverse datasets and/or species (e.g., rhesus to humans). We demonstrate the usefulness of transfer signatures in two use cases: the progression of latent to active tuberculosis and the severity of COVID-19 and influenza A H1N1 infection. This indicates that transfer signatures can be deployed in settings that lack disease-specific biomarkers. The broad significance of our work lies in the concept that a small set of archetypal human immunophenotypes, captured by transfer signatures, can explain a larger set of responses to diverse diseases.

Antibodies to the SARS-CoV-2 receptor-binding domain that maximize breadth and resistance to viral escape.

An ideal anti-SARS-CoV-2 antibody would resist viral escape 1-3 , have activity against diverse SARS-related coronaviruses 4-7 , and be highly protective through viral neutralization 8-11 and effector functions 12,13 . Understanding how these properties relate to each other and vary across epitopes would aid development of antibody therapeutics and guide vaccine design. Here, we comprehensively characterize escape, breadth, and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD), including S309 4 , the parental antibody of the late-stage clinical antibody VIR-7831. We observe a tradeoff between SARS-CoV-2 in vitro neutralization potency and breadth of binding across SARS-related coronaviruses. Nevertheless, we identify several neutralizing antibodies with exceptional breadth and resistance to escape, including a new antibody (S2H97) that binds with high affinity to all SARS-related coronavirus clades via a unique RBD epitope centered on residue E516. S2H97 and other escape-resistant antibodies have high binding affinity and target functionally constrained RBD residues. We find that antibodies targeting the ACE2 receptor binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency, but we identify one potent RBM antibody (S2E12) with breadth across sarbecoviruses closely related to SARS-CoV-2 and with a high barrier to viral escape. These data highlight functional diversity among antibodies targeting the RBD and identify epitopes and features to prioritize for antibody and vaccine development against the current and potential future pandemics.

Structural basis for broad sarbecovirus neutralization by a human monoclonal antibody.

The recent emergence of SARS-CoV-2 variants of concern (VOC) and the recurrent spillovers of coronaviruses in the human population highlight the need for broadly neutralizing antibodies that are not affected by the ongoing antigenic drift and that can prevent or treat future zoonotic infections. Here, we describe a human monoclonal antibody (mAb), designated S2X259, recognizing a highly conserved cryptic receptor-binding domain (RBD) epitope and cross-reacting with spikes from all sarbecovirus clades. S2X259 broadly neutralizes spike-mediated entry of SARS-CoV-2 including the B.1.1.7, B.1.351, P.1 and B.1.427/B.1.429 VOC, as well as a wide spectrum of human and zoonotic sarbecoviruses through inhibition of ACE2 binding to the RBD. Furthermore, deep-mutational scanning and in vitro escape selection experiments demonstrate that S2X259 possesses a remarkably high barrier to the emergence of resistance mutants. We show that prophylactic administration of S2X259 protects Syrian hamsters against challenges with the prototypic SARS-CoV-2 and the B.1.351 variant, suggesting this mAb is a promising candidate for the prevention and treatment of emergent VOC and zoonotic infections. Our data unveil a key antigenic site targeted by broadly-neutralizing antibodies and will guide the design of pan-sarbecovirus vaccines.

N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2.

The SARS-CoV-2 spike (S) glycoprotein contains an immunodominant receptor-binding domain (RBD) targeted by most neutralizing antibodies (Abs) in COVID-19 patient plasma. Little is known about neutralizing Abs binding to epitopes outside the RBD and their contribution to protection. Here, we describe 41 human monoclonal Abs (mAbs) derived from memory B cells, which recognize the SARS-CoV-2 S N-terminal domain (NTD) and show that a subset of them neutralize SARS-CoV-2 ultrapotently. We define an antigenic map of the SARS-CoV-2 NTD and identify a supersite (designated site i) recognized by all known NTD-specific neutralizing mAbs. These mAbs inhibit cell-to-cell fusion, activate effector functions, and protect Syrian hamsters from SARS-CoV-2 challenge, albeit selecting escape mutants in some animals. Indeed, several SARS-CoV-2 variants, including the B.1.1.7, B.1.351, and P.1 lineages, harbor frequent mutations within the NTD supersite, suggesting ongoing selective pressure and the importance of NTD-specific neutralizing mAbs for protective immunity and vaccine design.

Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity.

SARS-CoV-2 can mutate and evade immunity, with consequences for efficacy of emerging vaccines and antibody therapeutics. Here, we demonstrate that the immunodominant SARS-CoV-2 spike (S) receptor binding motif (RBM) is a highly variable region of S and provide epidemiological, clinical, and molecular characterization of a prevalent, sentinel RBM mutation, N439K. We demonstrate N439K S protein has enhanced binding affinity to the hACE2 receptor, and N439K viruses have similar in vitro replication fitness and cause infections with similar clinical outcomes as compared to wild type. We show the N439K mutation confers resistance against several neutralizing monoclonal antibodies, including one authorized for emergency use by the US Food and Drug Administration (FDA), and reduces the activity of some polyclonal sera from persons recovered from infection. Immune evasion mutations that maintain virulence and fitness such as N439K can emerge within SARS-CoV-2 S, highlighting the need for ongoing molecular surveillance to guide development and usage of vaccines and therapeutics.

N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2.

SARS-CoV-2 entry into host cells is orchestrated by the spike (S) glycoprotein that contains an immunodominant receptor-binding domain (RBD) targeted by the largest fraction of neutralizing antibodies (Abs) in COVID-19 patient plasma. Little is known about neutralizing Abs binding to epitopes outside the RBD and their contribution to protection. Here, we describe 41 human monoclonal Abs (mAbs) derived from memory B cells, which recognize the SARS-CoV-2 S N-terminal domain (NTD) and show that a subset of them neutralize SARS-CoV-2 ultrapotently. We define an antigenic map of the SARS-CoV-2 NTD and identify a supersite recognized by all known NTD-specific neutralizing mAbs. These mAbs inhibit cell-to-cell fusion, activate effector functions, and protect Syrian hamsters from SARS-CoV-2 challenge. SARS-CoV-2 variants, including the 501Y.V2 and B.1.1.7 lineages, harbor frequent mutations localized in the NTD supersite suggesting ongoing selective pressure and the importance of NTD-specific neutralizing mAbs to protective immunity.